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ABSTRACT 


Shipboard combat systems must be designed to withstand moderate to severe 
excitation induced by underwater explosion. Current specifications for combat 
systems shock qualifications are mandated in MIL-S-901D. Analyzing the differences 
and relationships between the predicted shock excitation, as derived from previous 
ship shock trials, and that shock excitation which is produced by the U.S. Navy 
Mediumweight Shock Machine required by MIL-S-901D, a proposed modification 
to the existing shock test procedure is presented which will better represent the 


shock phenomena experienced by combat systems exposed to underwater explosion. 
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I. INTRODUCTION 


Shipboard combat systems equipments must be designed to withstand severe 
shock excitations induced by underwater explosion, either conventional or nuclear 
in origin. The underwater explosion delivers violent forces to the ship in the form of 
an incident shock wave pressure, gas bubble oscillations, cavitation closure pulses and 
various reflection wave effects. These complex shock induced forces propagate 
through the ship to various combat systems equipments, severely damaging them 
unless they are designed and tested to withstand such violent excitations. The ability 
of a naval vessel to carry out its mission, after being subjected to an underwater 
explosion, depends on the survivability of these equipments. 

Current specifications for building ships and shipboard equipments contain the 
requirements for shock loading which must be met by the vendor of shipboard 
equipments. In general, all critical equipment is required to pass a series of shock 
tests which are outlined in Military Specification (MIL-S-901D), "Shock Tests, High 
Impact; Shipboard Machinery, Equipment and Systems, Requirements For." This 
document specifies the shock qualification test procedures which are required of all 
shipboard machinery, equipment and systems which must resist high impact 
mechanical shock. Three different shock test methods are outlined, these include 
shock testing by a lightweight shock machine, a mediumweight shock machine, or a 


floating platform barge. Selection of the shock test method depends on the item’s 


size and weight. All these devices deliver high impact mechanical shock excitation to 
items affixed to them. The purpose of these tests, again, is to determine the 
suitability and survivability of machinery and equipment for use during and after 
exposure to severe shock excitation which may occur in wartime. 

The response of combat systems equipments to underwater explosion 1s 
basically vibrational in nature. The equipment tends to vibrate at its fundamental 
natural frequency, or a low range of natural frequencies, when excited by the shock 
wave. The maximum amplitude of the vibration usually occurs after the shock wave 
passes the ship. The shock waveform is remarkably different at different levels within 
the ship, due largely to the ship’s structural and material characteristics which cause 
the shock waveform to lengthen in duration and decrease in frequency as it 
propagates upward through the ship. In essence, the ship acts as a low pass 
mechanical filter which alters the characteristics of the propagating shock wave from 
one possessing high frequency components to one that contains relatively low 
frequency components, as noted by Scavuzzo, Lam and Hill (1988). Figure 1 depicts 
the described phenomena. Thus, the study of shock qualification for combat systems 
equipments, which are usually located in upper levels of the ship, is a vibration 
problem in which relatively low frequency equipment support foundation excitations 
are observed. 

The U.S. Mediumweight Shock Testing Machine (MWSM), required by MIL-S- 
901D, is currently used for shock qualification of shipboard equipment ranging from 


about 250 to 6000 Ibs. This machine and its application is the primary focus of this 
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Figure 1. Acceleration Waveforms at Various Levels within the Ship. 
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study. The MWSM generates short-duration, high-impulse, high-frequency excitation 
which is transmitted directly to the mounted test object. This type of high frequency 
excitation waveform is significantly different from the actual waveforms that have 
been observed at various equipment locations during ship shock trials. The 
differences can be reduced by substituting a specially designed "tuned" test mounting 
fixture for the default mounting fixtures currently used to affix test items to the 
MWSM. A tuned mounting fixture, designed to respond at specific natural 
frequencies when excited by the MWSM, will provide a better simulation of the 
actual shock phenomena experienced by shipboard equipment. 

This study examines the. differences between the predicted shock waveform 
characteristics that can be observed in a ship shock trial, as reported for three 
representative pieces of equipment modeled in the Underwater Research Division 
of David Taylor Research Center (DTRC/UERD) DDG-51 Class Ship Pre-Shock 
Trial Analyses provided by Costanzo and Murray (1991), and those produced by the 
MWSM, as required by MIL-S-901D. The use of a "tuned" mounting fixture is 
proposed as a modification to existing shock test procedures. This proposal will 
afford a better representation of the actual shock phenomena experienced by surface 
shipboard combat systems equipments when they are qualified on the U.S. Navy 


Mediumweight Shock Machine. 


Il. BACKGROUND PRESENTATION 


A. U.S. NAVY HIGH-IMPACT SHOCK MACHINE FOR MEDIUMWEIGHT 
EQUIPMENT (MWSM) 


1. Development 

The need for shipboard equipment shock qualification was recognized 
during World War II when substantial damage to shipboard equipment resulted not 
from direct hit by a shell, but rather by the blast effects of explosions which occurred 
within the vicinity of the ship. The shock wave traveled through the structures within 
the ship causing excessive vibration and permanent deformation, which rendered vital 
combat equipment useless. 

In 1940, the first shock qualification test machine was developed by 
General Electric for the Navy. Called the Navy High-Impact Shock Machine for 
Lightweight Equipment, it was only papaUle of testing equipment which weighed up 
to 250 Ibs. The need for a machine to test heavier equipment was recognized and, 
in 1942, Westinghouse Electric Corporation developed the first Shock Machine for 
Mediumweight Equipment. It was capable of testing equipment which ranged from 
250 to about 4500 Ibs. Today, it remains virtually the same, however, the rating has 
been extended to handle equipments weighing up to 6000 Ibs and special equipment 


mounting fixtures are permitted, as outlined in MIL-S-901D. 


2. Description 

Perhaps the best description of the U.S. Navy High-Impact Shock Machine 
for Mediumweight Equipment (MWSM) can be found in the Naval Research 
Laboratory Report 7396 by Clements (1972). Paraphrasing his description, the 
MWSM is a hammer-anvil table apparatus, as noted in Figure 2. It consists of a 
hammer, weighing 3000 lIbs., which swings through an arc of up to 270 degrees. The 
hammer height is adjusted from a position 180 degrees away from the hammer 
impact and the total weight on the anvil table is used as an argument to determine 
this height. The hammer strikes the 4500 Ib. anvil table from below and imparts an 
upward, uniaxial acceleration and velocity to it. The anvil table has a 60 by 60 inch 
mounting surface upon which the test items are affixed. The entire anvil table 
apparatus is bolted to the machine’s foundation. These bolts permit the anvil table 
apparatus to travel up to 3 inches vertically after hammer impact. The table travel 
distance can be decreased by using pneumatic jacks to vertically reposition the anvil 
table. The machine is embedded in a massive concrete block resting on heavy coil 
springs which isolate it from its surroundings. The impacting surfaces of the hammer 
and anvil are fitted with spherical hardened-stee]l impact plates which render the 
collision elastic. 

MIL-S-901D mandates that test items will be mounted to the anvil table 
by a fixture in a manner characteristic of its designed shipboard orientation, along 
with any anticipated supporting structures which may mitigate the shock experience. 


The equipment and mounting fixture configuration cannot exceed 7400 Ibs. Steel 


channels are used to construct the mounting fixture for the test item. The number 
and type of channels to be used are specified in MIL-S-901D. Both equipment 
weight and distance between anvil table mounting holes determine the number of 
channels to be used for the mounting fixture. The specified configuration tends to 
keep the natural frequency of the test equipment-mounting fixture-anvil table system 
between 55 and 72 Hertz. This, as noted by Clements (1972), was not by design, but 
rather an effort to keep the maximum stress in the channels to less than 35,000 psi 


in a static acceleration field of 50 g’s. 


3. MWSM Shock Waveform 

The MWSM may be modeled quite simply as a mass-spring-damper system 
subjected to base excitation, as noted in Figure 3, and presented in Clements (1972). 
The system base excitation is provided by the hammer and anvil elastic impact which 
results in vertical motion of the anvil table. The equipment mounting fixture’s 
stiffness properties and the loss of Sea due to friction at bolted joints and 
imperfections in material design, substantiate this simplified model. More elaborate 
models may be required to describe and analyze intricate test structures, but 
meaningful results can be obtained with this model. 

The mechanical shock waveform afforded by the MWSM can be described 
by a velocity or acceleration waveform generated by the hammer and anvil impact, 
as noted in Clements (1972). Paraphrasing the description, the hammer and anvil 
impact produces a well defined half-sinusoidal acceleration pulse having an 


approximate duration of one millisecond. This pulse imparts an upward velocity and 





AVL TABLE 


Figure 2. The Navy High-Impact Shock Machine for Mediumweight 
Equipment (MWSM). The Dotted Line Shows Hammer Path. 
Courtesy Clements (1972). 


acceleration to the anvil table which continues until the table travel distance is 
achieved at the stops, some two to four milliseconds after impact. At this event, a 
new set of transients occur which may interfere with the motion established. This 
"table reversal" is followed by another transient which occurs when the anvil table 
comes to rest. In addition, the half-sinusoidal acceleration pulse excites a 750 Hertz 
longitudinal mode of the anvil table. This appears as a damped vibration that persists 
for about five cycles. Together, these events produce a very high-energy, high- 
frequency complex waveform with peak accelerations ranging from 220g to 580g, 
depending on the hammer height. 

Figure 4 depicts the peak anvil table accelerations versus hammer height 
and hammer impact velocity. The relationship between peak anvil table accelerations 
and associated hammer height and hammer impact velocity is linear, reinforcing the 
elastic impact argument. The "table reversal" acceleration pulse, occurring sometime 
later, will be somewhat smaller than the initial peak impact acceleration largely due 
to frictional factors. It follows that the transient acceleration pulse arising when the 
anvil table comes to rest will be even smaller than the "table reversal" transient 
acceleration pulse. 

Thus, the major features of the measured MWSM anvil table acceleration 
waveform can be described as a series of three half-sine acceleration pulses. The 
first, due to initial hammer-anvil impact, with a duration of one millisecond, followed 
by a second, oppositely directed, smaller pulse occurring sometime later at table 


reversal, then, lastly, an even smaller pulse when the anvil table comes to rest. The 


SIMPLIFIED MODEL OF MWSM 


EQUIPMENT 
MOUNTED TO 


MWSM 





EQUIPMENT 


ANVIL TABLE 


C= SYSTEM DAMPING 
K= SYSTEM STIFFNESS 


Figure 3. Simple Model MWSM. Courtesy Clements (1972). 
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latter two smaller pulses depend on the adjustable table travel distance. For analysis 
purpose, only the initial and, by consequence, severest, half-sine anvil table peak 
acceleration pulse will be simulated, for the complexities of structural damping 
influences and the effects table travel distances preclude the accurate simulation of 
the latter two smaller acceleration pulses. This first acceleration pulse will be used 


as the base excitation in analyzing the tuned fixture model response. 


B. DTRC/UERD DDG-51 CLASS DECK HOUSE PRE-SHOT TRIAL ANALYSES 

A transient shock analysis of the DDG-51 Class Deck House was conducted by 
the DIRC/UERD in erepatation for the forthcoming DDG-51 Class shock trial. The 
preliminary report, by Costanzo and Murray (1991), was obtained along with the 
predicted shock excitation histories and analyses for various weight combat 
equipments located on the 0-3 level of the DDG-51 Class Ship. This information was 
crucial in assessing the relationships between the shock phenomena experienced by 
surface shipboard equipment exposed to underwater explosion and that shock 


phenomena associated with the MWSM. Their findings are summarized below. 


1. Finite Element Model for DDG-51 Class Forward Deck House 
A finite element model of the DDG-51 Class Forward Deck House was 
developed which included all major structural members and supporting equipments. 
The model analysis was performed using COSMIC/NASTRAN code. Figure 5 
presents a depiction of the finite element model of the Forward Deck House. The 


model extends vertically from the 0-1 Level to the Sea Director Level. 
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Figure 4. MWSM Peak Anvil Table Accelerations Versus Hammer 
Height. Courtesy Clements (1972). 
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Figure 5. Finite Element Representation of DDG-51 Class Forward 


Deck House. Courtesy Costanzo and Murray (1991). 


The primary combat systems equipments present in the deck house were 
modeled by Costanzo and Murray (1991). The equipments were represented as 
lumped, rigid masses in the model. These masses were distributed to the various 
nodes present in the deck house model which represent the corresponding shipboard 
equipment locations in the actual deck house. In cases where the equipment center 
of mass was known, an appropriate offset was employed in the modeling procedure 
to position the equipment mass at the proper location above the respective deck or 
foundation attachments. Of particular interest were three combat systems equipments 
ranging from 325 to 4600 Ibs, all located on the 0-3 Level. The three included a 
Radar Receiver Transmitter (RT-1293/SPS-67) weighing 325 Ibs, a Beam 
Programmer (MX-10873\SPY-1D) richie 1000 Ibs and a Radio Frequency 
Amplifier (AM-7159\SP Y-1B) weighing 4600 lbs. Their range of weights would be 
useful in characterizing the behavior of the MWSM for various weight class 


equipments, low, medium, and high, when a tuned mounting fixture is applied. 


2. Analysis of DDG-51 Class Forward Deck House Model 
The finite element model was shock analyzed for maximum shock trial 
severity, shot four in a series of four underwater explosion shots. Transient shock 
response calculations for all nodes were performed using shock excitations to frames 
126, 174, and 220 at the 0-1 Level of the model, as noted in Figure 5. These three 
frames are the major supporting bulkheads of the DDG-51 Class Forward Deck 
House. The shock excitations were obtained from a full ship’s hull girder model of 


the DDG-51 class and actual shock trial data from previous cruiser shock trials. For 
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each significant deck house equipment, the dynamic responses at the equipment 
foundation locations were presented in the form of acceleration-time histories for the 
three principal directions, vertical, athwartships and fore/aft. For these directions, the 
transient shock responses for all three equipments were computed out to 70 
milliseconds. This was considered enough time for the ship to reach its maximum 
vertical displacement due to initial shock wave effect. 

Of primary importance to this thesis is an understanding of the Shock 
Spectra that results from these predicted transient acceleration excitations. The 
Shock Spectra defines the absolute maximum response envelope, over a wide range 
of system natural frequencies, of an undamped single degree of freedom mass-spring 
system subjected to a specific excitement. For a given excitement, the resultant Shock 
Spectra will reveal peak resonance responses which are of vital concern in the design 
and shock testing of equipment modeled as such a mass-spring system exposed to 
that excitement. 

Costanzo and Murray (1991) analyzed shock phenomena in three principal 
orientations: athwartships, fore/aft and vertical. Of particular interest is the vertical 
orientation shock analysis since this is the most severe type of shock experienced by 
surface vessels exposed to underwater explosion, as revealed in the DIRC/UERD 
Shock Spectra comparisons presented in Appendix A. Results of their vertical 
transient shock analysis, performed for the three equipments mentioned earlier, are 
shown in Figures 6 through 11. These are the predicted acceleration waveform 


excitements for each equipment’s foundation and the resultant Shock Spectras. It 
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should be noted that structural damping was omitted in their analysis for the 
following two reasons. One, damping is not constant throughout such a complex 
structure as the DDG-51 Class Forward Deck House and, two, the omission of 
structural damping generally results in conservative computed response levels. 

As noted in Figures 6 through 11, the predicted acceleration waveforms 
and associated Shock Spectras for the selected equipments, all located in the same 
compartment on the 0-3 Level, are dramatically different. Figures 6 through 8 shows 
the predicted equipment’s foundation acceleration excitement due to the shock wave. 
Each acceleration waveform is significantly different, possessing different frequency 
components and amplitudes. Placin g each of these equipment on the MWSM with 
the default mounting fixture outlined in MIL-S-901D and then exciting it with an 
acceleration pulse will not simulate the same shock phenomena depicted in those 
Figures. Thus, the need to apply a special MWSM mounting fixture, "tuned" to 
emulate the frequencies of interest, is necessary in order to provide the same 
characteristic shock phenomena observed in the field. 

The identification of the necessary characteristics of a tuned mounting 
fixture is the focus of this thesis. To this end, the three equipments studied provide 
an excellent representation of the weight ranges, low, medium and high, of the 
equipments tested on the MWSM. The findings presented in this thesis will enable 
the design, construction and implementation of tuned mounting fixtures on the 


MWSM. 
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Acceleration Time History for Radar Receiver/ Transmitter, 
RT-1293/SPS-67. Equipment Weight 325 lbs. Courtesy of 
Costanzo and Murray (1991). 


17 





z. 
& 
= 
< 
o 
= 2 
Z m 
a) a 
onl 
% 6 2 
el ¢ 
o = 
So & 
nile & 
~ Sf ” 
Pe — FE = 
o) 
3 © 
© 7 
Zz. 
oe {x3 
YW) 
Z = 
© = 
fm 
2 
=) 
fa) 
oe 
Ou. 
od 
= 
jon 
& 
fx) 
oa 
O, 
ve 
2) 
& 
a 
© © © © © © © 
ae) ~N = _ N 307 
| 
D 'NOILVYFTAIOV 
Figure 7. DDG-S51 Class Pre-Shock Trial Prediction of Vertical 


Acceleration Time History for Beam Programmer, MX- 
10873/SPY-1D. Equipment Weight 1000 lbs. Courtesy of 
Costanzo and Murray (1991). 
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Acceleration Time History for Radio Frequency Amplifier, 
AM-7158/SPY-1B. Equipment Weight 4600 Ibs. Courtesy of 
Costanzo and Murray (1991). 
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Spectra for Radar Receiver/Transmitter, RT-1293/SPS-67. 
Equipment Weight 325 Ibs. Courtesy of Costanzo and Murray 
(1991). 
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Figure 10. DDG-51 Class Pre-Shock Trial Predictions of Vertical Shock 


Spectra for Beam Programmer, MX-10873/SPY-1D. 
Equipment Weight 1000 Ibs. Courtesy of Costanzo and Murray 
(1991). 
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Figure 11. DDG-S51 Class Pre-Shock Trial Predictions of Vertical Shock 


Spectra for Radio Frequency Amplifier, AM-7159/SPY-1B. 
Equipment Weight 4600 lbs. Courtesy of Costanzo and Murray 
(1991). 
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Ill. SINGLE DEGREE OF FREEDOM TUNED MOUNTING FIXTURE 


A. DTRC/UERD PRE-SHOT TRIAL SHOCK SPECTRA ANALYSIS 

If the item to be tested is mounted to the MWSM by a special fixture, which 
has a designed fundamental frequency or frequencies when the item is affixed to it, 
it will experience the fixture’s vertical response to the MWSM excitation as its shock 
excitation, as presented by Chalmers and Shaw (1989). As the MWSM is an uniaxial 
machine, the tuned mounting fixture must be designed to provide a desired 
frequency response along a single direction, vertical for this study. The question then 
arises as to what fundamental frequency or frequencies to select. ‘The answer can be 
revealed by analyzing the item’s predicted Pre-Shot Trial Shock Spectra. 

The Vertical Orientation Shock Spectra, for each equipment provided by 
DTRC/UVERD, is presented in Figure 12. Each spectra shows a variety of peak 
responses at discrete frequencies. Beginning with the low weight range equipment 
Shock Spectra, a Radar Receiver/Transmitter weighing 325 Ibs, analysis shows two 
dominant peaks, one at about 60 Hertz and another at about 155 Hertz. The peak 
ratio is about 2:1. Clearly, this equipment’s predicted foundation acceleration 
excitement possesses two dominant waveform components. Thus, a two degree of 


freedom (DOF) uniaxial mounting fixture will be necessary in order to simulate the 
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dominant acceleration wavefonn characteristics present, as revealed within this Shock 
Spectra. 

Next, the medium weight range equipment Shock Spectra, a Beam Programmer 
weighing 1000 lbs, shows three well defined dominant peaks. One at about 23 Hertz 
and another at about 60 Hertz, both with about the same magnitude. An absolute 
dominant peak, by a factor of two, occurs at about 155 Hertz. A two DOF uniaxial 
tuned mounting fixture would be required as a minimum in order to simulate the two 
most important frequency characteristics depicted within this Shock Spectra. 

Finally, the high weight range equipment Shock Spectra, a Radio Frequency 
Amplifier weighing 4600 lbs, shows a single dominant peak at about 23 Hertz. A 
single DOF tuned mounting fixture would simulate the frequency and acceleration 


waveform characteristics found within this Shock Spectra. 


B. SINGLE DEGREE OF FREEDOM TUNED MOUNTING FIXTURE MODEL 

Advancing the principle proposed by Chalmers and Shaw (1989), the item to 
be tested is affixed to the MWSM by a single DOF mounting fixture, which has a 
designed fundamental frequency when the item is affixed to it. The test item will 
experience the fixture’s vertical response to the MWSM excitation as its shock 


excitation. Figure 13 depicts the described concept. 


1. Single Degree of Freedom Tuned Mounting Fixture Mathematical Model 
Figure 13 depicts a single DOF mass-damper-spring system subjected to 


foundation excitement, the model representing the single DOF tuned mounting 
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DTRC/UERD Predicted DDG-51 Class Pre-Shock Trial 


Figure 12. 
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Vertical Orien 
Murray (1991). 


fixture. The equation of motion for this system can be expressed in terms of relative 
motion coordinates. This will facilitate the solution to the problem. The mass in this 
model represents the combined equipment and mounting fixture mass. Let the 
absolute motion of the mass He expressed by the x coordinate and the foundation 


motion by the z coordinate. Then equations 1 through 3: 


y=x-z (1) 
y=t-2 (2) 
p=k-Z (3) 


are the relative coordinate transformations for displacement, velocity and 
acceleration, respectively. 

The system’s natural frequency is a function of system mass, m, and 
stiffness, k, and can be expressed either in radians per second or in Hertz, cycles per 


second. Equations 4 and 5 note those respective relationships. 
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SINGLE DEGREE OF FREEDOM 
TUNED MOUNTING FIXTURE 
MODEL 





M = SYSYEM MASS — 
K = SYSTEM STIFFNESS 
C = SYSTEM DAMPING 


Figure 13. Single Degree of Freedom Tuned Mounting Fixture. 
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The equation of motion for this system can be written as follows: 


J+2Ew yrwry=-Z (6) 


where € is the damping factor and must be estimated. 

Equation 6 was numerically integrated using an unconditionally stable 
numerical integration scheme outlined by Craig (1981). The Fortran code listing is 
presented in Appendix B. 

Once the relative displacements, velocities, and accelerations are known, 


equations 1 through 3 can be solved for the absolute quantities: 


X=y+z (7) 
X=Y+Z (8) 
X=f+7Z (9) 
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where equation 9 presents the single degree of freedom tuned mounting fixture’s 
response acceleration, which is the equipment’s foundation excitement. This 
acceleration excitement is used for the development of the Shock Spectra. The 


Fortran code listing for the Shock Spectra formulation is listed in Appendix C. 


2. Single Degree Of Freedom Tuned Mounting Fixture Modeled Application 

The high weight range equipment, a Radio Frequency Amplifier weighing 
about 4600 Ibs, was an excellent candidate for the single DOF tuned mounting 
fixture application. The DIRC/UERD Shock Spectra revealed a peak at about 23 
Hertz, thus, the mounting fixture must possess a fundamental frequency of 23 Hertz 
when the equipment is affixed to it. This can be easily obtained by first, fixing the 
mass of the system, then designing the system stiffness, noting the relationships 
presented in equations 4 and 5. 

As stated earlier, the major feature of the MWSM excitation 1s a half-sine 
acceleration pulse of approximately one millisecond duration. Selection of MWSM 
hammer height and consequent peak anvil table acceleration is dependent on total 
anvil table top weight, as required by MIL-S-901D. For the single DOF tuned 
mounting fixture application analysis, a peak acceleration of 425g’s was selected. This 
corresponds to a hammer height of three feet, which is nominal for the first series 
of hammer blows corresponding to total anvil table top weights in the range of 4600 
to 7400 Ibs. Figure 14 depicts the simulated MWSM pulse used for the single DOF 


tuned fixture application analysis. 
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A sensitivity study in damping was conducted to correlate and validate the 
single DOF tuned mounting fixture’s modeled response with respect to actual 
MWSM calibration test data. Examination of MWSM calibration data, as compiled 
by Costanzo and Clements (1988), showed that the test weight acceleration responses 
appeared to dissipate within .2 to .4 seconds after initial MWSM excitement. Thus, 
the selection of the damping factor in equation 6 could be estimated based on that 
information. Figures 15 through 19 show the damped acceleration response study for 
the single DOF tuned mounting fixture subjected to the MWSM acceleration pulse 
excitation. 

Figure 18, with a damping factor equal to .08, presents the best decaying 
characteristics consistent with the MWSM test calibration data. The single DOF 
tuned mounting fixture, damped at this factor, was further processed to yield the 
acceleration waveform and comparative results presented in Figures 20 to 24. Figure 
20, the Fourier Transform of the acceleration response, shows a well defined peak 
at 23 Hertz, as expected. Figure 21 shows the comparison between the DIRC/UERD 
predicted acceleration waveform excitement for the Radio Frequency Amplifier and 
that produced in the first 70 milliseconds by the single DOF tuned mounting fixture. 
There is close agreement in shape and magnitude, as expected. 

Figure 22 is the resultant Shock Spectra, which is typical for a decaying 
sinusoidal acceleration excitement. The modeled single DOF tuned mounting fixture 
Shock Spectra is compared with the predicted DTRC/UERD Shock Spectra in 


Figure 23. There is excellent agreement in spectral shape. The magnitude difference 
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SDOF TUNED FIXTURE RESPONSE TO MWSM ACCELERATION PULSE 
SECONDS 





Figure 15. Modeled Single DOF Damped Tuned Mounting Fixture 


Acceleration Response to MWSM Acceleration Pulse of Peak 
425g’s. 
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SDOF TUNED FIXTURE RESPONSE TO MWSM ACCELERATION PULSE 
SECONDS 





Figure 16. Modeled Single DOF Damped Tuned Mounting Fixture 
Acceleration Response to MWSM Acceleration Pulse of Peak 
425¢’s. 
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SDOF TUNED FIXTURE RESPONSE TO MWSM ACCELERATION PULSE 
SECONDS 





Figure 17. Modeled Single DOF Damped Tuned Mounting Fixture 


Acceleration Response to MWSM Acceleration Pulse of Peak 
4258's. 
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SDOF TUNED FIXTURE RESPONSE TO MWSM ACCELERATION PULSE 
SECONDS 





Figure 18. Modeled Single DOF Damped Tuned Mounting Fixture 
Acceleration Response to MWSM Acceleration Pulse of Peak 
4259’s. 
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SDOF TUNED FIXTURE RESPONSE TO MWSM ACCELERATION PULSE 
SECONDS 





Figure 19. Modeled Single DOF Damped Tuned Mounting Fixture 
Acceleration Response to MWSM Acceleration Pulse of Peak 
425¢’s. 
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arises due to the duration of excitement. For this study, the complete decaying 
sinusoidal] base excitement was allowed in the formulation of the Shock Spectra. In 
the DTRC/UERD study, only 70 milliseconds of undamped base excitation was 
permitted, as this would allow the ship to achieve maximum vertical displacement 
due to initial shock excitation and, as stated earlier, damping influences were 
omitted. Closer agreement in shock spectra shape and magnitude at resonance would 
result with a longer duration of base excitement in the pre-shock trial analysis, after 
the complexities of damping are investigated. Figure 24 shows a the response of an 
equipment, modeled as an undamaged single DOF system, subjected to the tuned 
fixture’s acceleration excitation. 

The development of an MWSM single DOF tuned mounting fixture has 
been under investigation for sometime, as noted by Chalmers and Shaw (1989). The 
Soft Deck Simulator, shown in Figure 25, is such a device. The Soft Deck Simulator, 
developed by Naval Underwater Systems Center, was intended for the shock 
qualification of submarine combat systems equipments on the MWSM. It is a simple 
single DOF mass-spring system composed of springs, which are inserted in a parallel 
configuration by cartridges. Each cartridge possesses 12 springs which are sandwiched 
between support rails. The total weight of springs, support rails and test equipment 
must be considered when selecting the frequency response characteristics of the 
system. By varying the number of cartridges, the desired frequency response can be 


achieved. It is capable of achieving frequency responses in the range of about 19 to 
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30 Hertz. By the above analysis, this device could be implemented for the shock 


qualification of surface combat systems equipments with very promising results. 
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FIXTURE TUNED TO 23 HZ 
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FFT SDOF TUNED FIXTURE ACCELERATION RESPONSE 
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Figure 20. Fourier Transform of Single DOF Damped Tuned Mounting 
Fixture Acceleration Response. (Zeta = .08) 
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Pre-Shock Trial for Node 3320, Weight 4600 Ibs, and Single 


DOF Tuned Mounting Fixture. 


Figure 21. 
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Comparative Shock Spectra. Single DOF Damped Tuned 
42 


Mounting Fixture and DIRC/UERD Predicted Pre-Shock 


Tnial for Frequency Amplifier Weight 4600 Ibs. 


Figure 23. 
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IV. TWO DEGREE OF FREEDOM TUNED MOUNTING FIXTURE 


A. DTRC/UERD PRE-SHOT TRIAL SHOCK SPECTRAL ANALYSIS 

As presented earlier, both the light weight range and medium weight range 
equipment possess Vertical Orientation Shock Spectra that reveal two or more well 
defined dominant peaks at discrete frequencies. The light weight range equipment, 
a Radar Receiver/Transmitter weighing 325 Ibs, possesses two peaks, one at 60 and 
the other at 155 Hertz, with a 2:1 magnitude ratio. A two DOF uniaxial tuned 
mounting fixture would provide the frequency characteristics necessary to simulate 
this phenomena. 

The medium weight range equipment, a Beam Programmer weighing 1000 Ibs, 
possesses three well defined dominant peaks. Two peaks, with about the same 
Magnitude, occur at about 23 and 60 Hertz, respectively. The third and absolute 
dominant peak, by a factor of two, occurs at about 155 Hertz. A two DOF uniaxial 
tuned mounting fixture could simulate the two most important frequency 
characteristics revealed within this Shock Spectra. Knowing that higher accelerations 
are experienced by equipment possessing higher fundamental frequencies, the 60 
Hertz and 155 Hertz frequencies are deemed the two most important frequency 
characteristics for this study. Further investigation of other combinations are 


possible. 
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B. TWO DEGREE OF FREEDOM TUNED MOUNTING FIXTURE MODEL 
Extending the concept developed by Chalmers and Shaw (1989) to a two DOF 
tuned mounting fixture is the focus of this chapter. Figure 26 depicts the extended 
concept. The upper tier of the system possesses the test item of interest. The coupled 
response of this tier to base shock excitation will be the shock excitement 
experienced by the item tested. The mass, damping and stiffness relations of the two 


DOF tuned mounting fixture model will need to be investigated. 


1. Two Degree of Freedom Tuned Mounting Fixture Mathematical Model 
Referring to figure 26 again, let the upper tier, with mass m,, be comprised 
of the equipment tested along with its associated support mountings. The lower tier, 
with mass m,, is comprised of test weights and support mountings. Each tier 
possesses characteristic damping and stiffness properties, c and k, respectively. 
Expressing the absolute coordinates, x, and x,, in terms of relative coordinates yields: 
y, =X, -Z 


y,=%,-2 (10) 
Y, =, -2 


=X -2 


yak <8 (11) 
Yq =X_~Z 


where the z coordinate represents the foundation motion. Once again, if the relative 
coordinate quantities are known, then it is a simple matter to obtain the absolute 


quantities. 
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Figure 26. Two Degree of Freedom Tuned Mounting Fixture Model. 
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The coupled equations of motion for this system can be obtained and are 


expressed in matrix form below. 
yy E my 0 =Z (12) 
5 O m,||-z 


k, -k, 


“ky (ky +k) 





m, ONY, C) “C, Wy 
ee a e 7 
0 m,))y¥, “€, (€,+C,) ||¥, 


Each tier has its own natural frequency which may be expressed in either 











radians per second or in Hertz, cycles per second. Equations 13 through 16 present 
those relationships. 


£ = (13) 
fot = (14) 
a : (15) 

fo 2 (16) 


If one were to separate the tiers, each could be viewed as a single DOF 


mass-damper-spring system and the overall system response could be formulated by 
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modal analysis. Equations 14 and 16 are the uncoupled natural frequencies of each 
tier. For this study, the coupled natural frequencies are of interest, since it is at those 
frequencies that the system response is observed. Defining the ratio of the upper tier 


mass to the lower tier mass as follows: 


mass, (17) 





Mass, 


it can be shown, as noted in Shin (1981), that the relationship between the system’s 


coupled natural frequencies and uncoupled frequencies are: 


fry 8h fe AG ah OLE PV (18) 


Fras oh Uh 8hs sh) Of FV (19) 


where fn, and fn, are the system’s coupled natural frequencies in Hertz, which are 


easily converted to radians per second as noted below. 


on, =2nfn, (20) 
Wn, =2nfn, (21) 
Knowing both the system’s coupled and uncoupled natural frequencies, 


along with the system’s mass ratio, equation 12 may be decoupled by modal analysis, 


49 


as presented in Mierovitch (1986), and solved as two single DOF systems. Once the 
response of each is single DOF is known, the coupled response may be obtained. 
Equations 22 and 23 show the decoupled equations of motion using the coupled 
natural frequencies, wn, and wn,. It must be noted that the system damping is 
assumed to be a linear combination of system mass and stiffness. This valid 


assumption will permit the modal analysis approach. 
q) +2& omg, +n, q,=-%, (22) 
G,+28 wn,g,+wn,q,=-2, Es) 


The natural coordinates, q, and q,, are related to the relative coordinates, 


y, and y,, by the modal matrix as follows: 
uy, u 
p a 1] 7 i (24) 
ooi| bbe 


Equations 22 and 23 were solved with an unconditionally stable numerical 
integration scheme and then the absolute motions were obtained per the relations 
in equations 10 and 11. The Fortran code listing for this numerical solution is 


presented in Appendix D. 
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2. Two Degree of Freedom Tuned Mounting Fixture Modeled Application 

The DTRC/UERD shock spectral analysis revealed two common 
frequencies of interest between the light weight and medium weight range 
equipment. Those frequencies were 60 and 155 Hertz. These are the coupled 
response frequencies that must be observed in the modeled application. As stated 
earlier, knowing the two uncoupled system natural frequencies, f, and f,, along the 
mass ratio, the coupled natural frequencies, fn, and fn,, may be obtained by 
equations 18 and 19. However, only the coupled frequencies, fn, and fn, are known, 
thus, a sensitivity study is required in order to attain reasonable values for the 
uncoupled tier frequencies, f, and f,, and system mass ratio. 

Figures 27 through 31 present the results of a study performed where the 
lower tier natural frequency, f,, was held constant and iterations of f, were conducted 
to ascertain the combination of f,, f, and mass ratio that would provide the desired 
coupled frequency response at 60 and 155 Hertz. Mass ratios of 1 and .1 were 
studied. In all cases, the coupled frequencies for a mass ratio of .1 were bounded by 
those of mass ratio 1, thus defining a mass ratio coupled frequency envelope. Figure 
30 shows that, at a mass ratio of 1 and f, equal to 100 Hertz, a selection of f, equal 
to 94 Hertz would provide the desired coupled frequency response for this study. In 
addition, a mass ratio of 1 would minimize the total anvil table top weight, thus 
permitting a wider application of the two DOF tuned mounting fixture. For the 


remainder of this analysis, the upper tier natural frequency is 94 Hertz and the lower 


a1 


tier possesses a natural frequency of 100 Hertz. The Fortran code for this iterative 
scheme is provided in Appendix E. 

Knowing both the set of coupled and uncoupled natural frequencies, along 
with a mass ratio of 1, an analysis of the two DOF tuned mounting fixture was then 
permitted. Tier weights of 1700 lbs each were assumed. This would allow a wide 
range of possibilities for equipment and support mounting combinations and would 
result in a moderate total anvil table top weight of 3400 lbs. Further, knowing the 
respective tier mass and natural frequency properties, the equivalent tier stiffnesses 
could be calculated. 

The MWSM acceleration half-sine acceleration pulse for this analysis is 
presented in Figure 32. As this model is a relatively stiff system where higher 
accelerations can be expected, the peak MWSM acceleration pulse had to be 
adjusted to provide meaningful results. A peak acceleration of 75g’s was used in the 
modeled application analysis of the two DOF tuned mounting fixture. The MWSM 
is capable of delivering such a peak acceleration. A MWSM calibration study, based 
on the analysis of the two DOF tuned mounting fixture, will be necessary to 
determine suitable hammer heights for shock testing equipment mounted to the 
machine by this fixture. 

The same type of sensitivity study in damping was conducted as in the 
single DOF tuned mounting fixture analysis. For simplicity, both tiers were subjected 
to the same damping factor, however, other combinations are possible. Figures 33 


through 36 show the effects of damping on the upper tier acceleration response. 
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Figure 27. 
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Figure 28. 
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Figure 29. 
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Figure 30. 
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Figure 31. 
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Figure 34 , with a damping factor of .02, shows the best decaying characteristics 
consistent with the MWSM calibration data, compiled by Costanzo and Clements 
(1988), and was used in the further analysis of the two DOF tuned mounting fixture. 

Figure 37 shows the Fourier Transform of the upper tier acceleration. As 
expected, there are two peaks, one at 60 Hertz and the other at 155 Hertz. The 
magnitude of each peak is proportional to the amplitude of that frequency 
component in the acceleration waveform. The dominant first frequency component 
indicates that the first mode of vibration has a greater contnbution to the response. 
Figures 38 and 39 show the comparisons of the first 70 milliseconds of the modeled 
upper tier’s acceleration waveform with those predicted by the DIRC/UERD study 
for the low weight and medium weight range equipment. Very close agreement exists 
between the modeled acceleration waveform and the acceleration waveform for the 
low weight range equipment. Further studies as to the effects of tier weights on the 
acceleration waveforms will provide closer agreement between the modeled upper 
tier waveform and that waveform predicted by a pre-shock trial] analysis. 

Figure 40 shows the resultant Shock Spectra using the upper tier’s 
acceleration response as base excitation. The shape and relative peak magnitudes are 
typical for such an excitement. Figures 41 and 42 show the comparison between the 
DTRC/UERD predicted Shock Spectra and that resultant from this study of a two 
DOF tuned mounting fixture. There is very close agreement in apectral shape and 
magnitude with respect to the low weight range equipment’s Shock Spectra. Closer 


agreement in the Shock Spectra of the medium weight range equipment will result 
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after further studies are conducted as to the effects of tier weights in the coupled 
response. Further, closer agreement would result in the spectras if longer duration 
base excitation were permitted in the pre-shock trial analysis. This, as noted 
previously, can occur once the complexities of damping are investigated, as noted 
previously. Figure 43 shows the typical resonance response of an equipment, modeled 
as an undamped single DOF system, subjected to the upper tier’s decaying sinusoidal 
acceleration, a waveform consisting of two frequency components, one at 60 Hertz 
and the other at 155 Hertz. 

Development of a two DOF tuned mounting fixture with a coupled 
frequency response of 60 and 155 Hertz can be obtained once the mass ratio is 
selected. For this study, a mass ratio of 1 optimized the total anvil table top weight 
and provided reasonable tier natural frequencies of 94 and 100 Hertz. To achieve 
those relatively high tier natural frequencies, the proposed model in Appendix F can 
be constructed. The equivalent stiffness properties of each tier are provided by the 
tier support mounting beam configurations and the manner in which those beams are 
loaded. A study into the design and construction of such a mounting fixture is 


warranted by the above analysis. 
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Figure 32. Simulated MWSM_ Half-Sine Acceleration Pulse for 


Evaluation of the Two DOF Tuned Mounting Fixture. 


2 DOF TUNED FIXTURE UPPER TIER ACCELERATION RESPONSE 
SECONDS 





Figure 33. Two DOF Tuned Mounting Fixture Upper Tier Damped 
Acceleration Waveform. Zeta = .01 
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2 DOF TUNED FIXTURE UPPER TIER ACCELERATION RESPONSE 
SECONDS 





Figure 34. Two DOF Tuned Mounting Fixture Upper Tier Damped 
Acceleration Response. Zeta = .02. 
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Figure 35, Two DOF Tuned Mounting Fixture Upper Tier Damped 


Acceleration Response. Zeta = .03. 
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2 DOF TUNED FIXTURE UPPER TIER ACCELERATION RESPONSE 
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Figure 36. Two DOF Tuned Mounting Fixture Upper Tier Damped 


Acceleration Response. Zeta = .04. 
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Figure 38. DTRC/UERD Predicted Acceleration Waveform for Node 


3310, Radar Receiver/Transmitter, and Two DOF Tuned 
Mounting Fixture Upper Tier Acceleration Waveform for First 
70 Milliseconds. 
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Figure 39. 


DTRC/UERD PRETRIAL PREDICTION, NODE 3314 VERTICAL ACCELERATION 
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DTRC/UERD Predicted Acceleration Waveform for Node 
3314, Beam Programmer, and Two DOF Tuned Mounting 
Fixture Upper Tier Acceleration Waveform for First 70 
Milliseconds. 
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Figure 40. Shock Spectra Using Two DOF Tuned Mounting Fixture 
Upper Tier Acceleration as Excitation. 
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SHOCK SPECTRA: 2 DOF TUNED FIXTURE VS DTRC/UERD PREDICTED 


Figure 41. 
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Shock Spectra Comparison. DTRC/UERD Predicted for Node 
3310, Radar Receiver/Transmitter, Foundation Excitation and 
Upper Tier of Two DOF Tuned Mounting Fixture Excitation. 
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Figure 42. 


SHOCK SPECTRA: 2 DOF TUNED FIXTURE VS DTRC/UERD PREDICTED 
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Shock Spectra Comparison. DIRC/UERD Predicted for Node 
3314, Beam Programmer, Foundation Excitation and Upper 
Tier of Two DOF Tuned Mounting Fixture Excitation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The DIRC/UERD DDG-51 Class Pre-Shock Trial Shock Analyses of three 
combat systems equipments ranging from 325 to 4600 Ibs revealed shock 
characteristics that cannot be simulated on the MWSM unless special tuned 
mounting fixtures are implemented. Both the single DOF and two DOF tuned 
mounting fixtures can be used to produce a more realistic shock phenomena when 
test items are shock qualified under the simulated conditions inherent with the 
MWSM. The proposals below are presented for consideration in advancing the 
practice of shock qualifying surface ship combat systems equipments. 

For heavy weight equipment, in the range of about 4600 Ibs, primarily low 
frequency foundation excitation can be expected to dominate, yielding acceleration 
waveforms consisting of one dominant frequency component in the vicinity of about 
23 Hertz. A single DOF tuned mounting fixture used to affix a test item to the 
MWSM will provide the shock characteristics observed in this situation. The Soft 
Deck Simulator, developed by the Naval Underwater Systems Center for the shock 
qualification of submarine combat systems equipments, is proposed for use in the 
shock qualification of heavy weight range surface ship combat systems equipments 
which display the foundation excitations described above. Analysis of the ship class 
pre-shock trial data will reveal which equipments are likely to experience support 


foundation excitations that can be simulated by such a device. 


g2 


For low and medium weight equipments, in the ranges of 325 and 1000 lbs, 
respectively, the DIRC/UERD analyses revealed a more complex waveform in which 
a two DOF tuned mounting fixture must be used. It is proposed that the two DOF 
tuned mounting fixture, described in chapter IV, be used to simulate the dominant 
shock characteristics revealed in the ship class pre-shock trial analyses for low and 
medium weight range combat systems equipments. It is strongly recommended that 
the below proposals be implemented to advance the practice of shock qualifying 
surface ship combat systems equipments: | 

@ the design and construction of a two DOF tuned mounting fixture for the 

MWSM 

® a sensitivity/calibration study of tier weights and frequency response 
characteristics of the two DOF tuned mounting fixture for given MWSM 


hammer heights 


® an investigation of damping characteristics in shock wave propagation through 
shipboard structures to advance the study of pre-shock trial analyses 


® the development of a weight category, low and medium, combat system 
equipment Shock Spectra library, for various vessel classes, as a reference in 
parameter selection for the two DOF tuned mounting fixture 
Implementation of the above proposals will produce a MWSM two DOF tuned 
mounting fixture with the ancillary information necessary to provide a more realistic 
shock phenomena when low and medium weight range combat systems equipments 
are shock qualified in a simulated environment, such as with the MWSM. 


This overall study demonstrated that the use of tuned mounting fixtures on the 


MWSM can be used to accurately simulate the shock characteristics that may be 
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observed in actual ship shock trials. Implementing these devices in the U.S. Navy 
Shock Qualification Program for Surface Ships will promote system reliability in 


times that are crucial to vessel survivability--the time of an underwater attack. 
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APPENDIX A 


The DTRC/UERD DDG-51 Class Pre-Shock Trial Analyses Shock Spectra for 
three modeled combat systems equipment located on the 0-3 Level: 


TABLE A-1 


| EQUIPMENT FINITE ELEMENT EQUIPMENT WEIGHT 
NODE LBS 
| RADAR RECEIVER/ 3310 
iardiea iene — 


BEAM | BEAM PROGRAMMER ~ Lr ti(‘(a<i‘CéCrS*? 
RADIO FREQUENCY 3320 fe 10 
AMPLIFIER 


Figures A-1 through A-3 reveal the three Shock Spectra orientations, fore/aft, 














athwartships and vertical for each equipment listed above. By comparison, the 
vertical orientation presents the severest shock phenomena experienced by each 


equipment. 
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DTRC/UERD Radar Receiver/Transmitter Shock Spectra. 
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Courtesy of Costanzo and Murray (1991). 


Figure A-1. 
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DTRC/UERD Beam Programmer Shock Spectra. Courtesy of 
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Costanzo and Murray (1991). 


Figure A-2. 
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DTRC/VERD Radio Frequency Amplifier Shock Spectra. 


Courtesy of Costanzo and Murray (1991). 


Figure A-3. 
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APPENDIX B 


C LT RANDALL CORBELL 

C NAVAL POSTGRADUATE SCHOOL 

C SHOCK QUALIFICATION OF COMBAT SYSTEMS EQUIPMENTS USING TUNED 

C MOUNTING FIXTURES ON THE U.S. NAVY MEDIUMWEIGHT SHOCK MACHINE 


C PROGRAM: SDOF TUNED FIXTURE RESPONSE TO MWSM HALF-SINE 
ec ACCELERATION PULSE 


: 


(A) CLEMENTS,E.W.,"SHIPBOARD SHOCK AND NAVY DEVICES FOR 
ITS SIMULATION",NRL REPORT 7396, 1972. 

(B) CRAIG,R,"STRUCTURAL DYNAMICS, AN INTRODUCTION TO 
COMPUTER METHODS", JOHN WILEY AND SONS, 1981. 

(C) COSTANZO/MURRAY, "DTRC/UERD DDG-51 CLASS PRESHOCK 
TRIAL ANALYSIS PRELIMINARY REPORT’,1991. 

(D) COSTANZO/CLEMENTS,"MWSM CALIBRATION DATA",1988. 


THIS PROGRAM COMPUTES THE RESPONSE OF A SINGLE DEGREE OF 
FREEDOM MASS-DAMPER-SPRING SYSTEM TO A BASE EXCITEMENT HALF-SINE 
ACCELERATION PULSE OF 1IMSEC DURATION, AS CHARACTERISTIC WITH THE 
MWSM NOTED IN REF(A). REF(B) WAS USED FOR NUMERICAL ANALYSIS 
PROGRAM STRUCTURE. THE ACCELERATION RESPONSE OF THE SDOF TUNED 
FIXTURE IS THE BASE EXCITATION A WHICH COMBAT SYSTEM EQUIPMENT 
WILL EXPERIENCE IN UNDEX, AS SIMULATED ON THE MWSM. 


SELECTION OF THE NATURAL FREQUENCY OF THE SYSTEM IS 

BASED ON ANALYSIS OF FFT/SHOCK SPECTRAL INFORMATION PROVIDED IN 
REF (C). SELECTION OF DAMPING VALUES WERE OBTAINED AFTER ANALYSIS 
OF MWSM CALIBRATION DATA, NOTED IN REF (D). 


DECLARATIONS... 
PARAMETER(MAX= 10000) 
ARRAYS... 
REAL A(3,3),L(3),UO(3),UI(3),U(3) 
DIMENSION TIME (MAX),ZBASE (MAX), TUNFIXACC(MAX), 
RELDIS(MAX),RELVEL(MAX),RELACC (MAX) 
VARIABLES... 
REAL ZETA,DELT,FREQ,E,F,G,WN,WN2,Z 
C INTEGERS... 
INTEGER NSTEP 


CO SEO Oa Ga a Ore) 6) Cy). 07) 


©) 


Cc PRINT STATEMENTS FOR INPUT... 
PRINT*, INPUT VISCOUS DAMPING COEFFICIENT,ZETA.’ 
READ*, ZETA 

Cc PARAMETERS USED... 
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FREQ = _ 23. 
DELT = 0.0001 
NUMSTEP = 10000 
UO(1l) = 0. 
UO(2) =O. 
REWIND 10 


OPEN(10,FILE=’ANVILACCHV.DAT ,STATUS=’OLD’) 
OPEN(20,FILE=’SDOFTUNFIX.DAT ,STATUS=’NEW’) 


c READ DATA FROM ANVIL ACCELERATION RECORD, NOTE ACCELERATION IN 


DO 100 I=1,NUMSTEP 
READ(10,*) TIME(I),ZBASE(I) 
100 CONTINUE 


C | SOLVE EQN MOTION FOR INITIAL ACCELERATION: KNOW INITIAL 
C  MASS,ZETA,FREQ,CONDITIONS... 
WN = 2.*4.*ATAN(1.)*FREQ 
WN2 = =(2.*4.*ATAN(1.)*FREQ)**2 
UO(3) =-ZBASE(1)*32.2-(2.*ZETA*WN)*UO(2)-WN2*UO(1 
Ul(1) =UO(1) | 
Ul(2) =U0(2) 
Ul(3) =U0(3) 
RELDIS(1) =UO(1) 
RELVEL(1) =UO(2) 
RELACC(1) =UO(3) 


€ COMPUTE PARAMETERS FOR OPERATOR MATRICES... 
E=WN2*(DELT**2) 
F=2.*ZETA*WN*DELT 
G=(1./(1.+(F/2.)+(E/4.))) 


C LOAD AMPLIFICATION AND LOAD VECTOR MATRIX... 
A(1,1)=G*(1.+(F/.)) 
A(1,2)=G*DELT*(1.+ (F/4.)) 
A(1,3)=G*DELT**2/4. 
A(2,1)=G*(-E/(2.*DELT)) 
A(2,2)=G*(1.-E/4.) 
A(2,3)=G*DELT/2. 
A(3,1)=G*(-E/((DELT* *2)) 
A(3,2)=G*(-1.)*(E+F)(DELT) 
A(3,3)=G*(-1.)*((F/2.) + (E/4.)) 
L(1) =G*(DELT**2)/4. 

L(2) =G*DELT/2. 
L3) =G 


C BEGIN ITERATIONS... 
DO 200 J=1,NUMSTEP 
Z =-ZBASE(J+1)*32.2 
U(1)= A(1,1)*UI(1)+A(1,2)*UI(2)+A(1,3)*U1(3)+Z*L(1) 
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22 
200 


250 


U(2)= A(2,1)*UI(1)+A(2,2)*U1(2) +A(2,3)*U1(3)+Z*L(2) 
U(3)= A(3,1)*UI(1)+A(3,2)*UI(2) +A (3,3)*UI(3) +Z*L(3) 


SAVE THE DIS, VELACC ARRAY VALUES FOR THIS ITERATION... 
RELDIS(J+1) = U(1) 
RELVEL(J+1) = U(2) 
RELACC(J+1) = U(3) 


RESET THE ITERATION MATRIX FOR THE NEXT TIME 

ITERATION... 

DO 225 K=1,3 
UI(K)=U(K) 

CONTINUE 

CONTINUE 

DO 250 I=1,NUMSTEP 
TUNFIXACC(I)=(RELACC(I)/32.2) +ZBASE(I) 
WRITE(20,*)TIME(I), TUNFIXACC(I) 

CONTINUE 

END 
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APPENDIX C 


LT RANDALL CORBELL 

NAVAL POSTGRADUATE SCHOOL 

SHOCK QUALIFICATION OF COMBAT SYSTEMS EQUIPMENTS USING TUNED 
MOUNTING FIXTURES ON THE U.S. NAVY MEDIUMWEIGHT SHOCK MACHINE 


Cl wae.) @) 


PROGRAM: SHOCK SPECTRA FORMULATION USING TUNED MOUNTING 
IXTURE 


29 


ACCELERATION EXCITATION 


REF: (A) CLEMENTS,E.W.,"SHIPBOARD SHOCK AND NAVY DEVICES FOR 
ITS SIMULATION",NRL REPORT 7396, 1972. 
(B) CRAIG,R,"STRUCTURAL DYNAMICS, AN INTRODUCTION TO 
COMPUTER METHODS", JOHN WILEY AND SONS, 1981. 
(C) COSTANZO/MURRAY, "DTRC/UERD DDG-51 CLASS PRESHOCK 
TRIAL ANALYSES PRELIMINARY REPORT", 1991. 


THIS PROGRAM COMPUTES SHOCK SPECTRA FOR AN UNDAMPED A 

SINGLE DEGREE OF FREEDOM MASS-SPRING SYSTEM SUBJECTED TO A BASE 

EXCITEMENT, THE TUNED MOUNTING FIXTURE ACCELERATION ON THE 
WSM 

IN REF(A). REF(B) WAS USED FOR NUMERICAL ANALYSIS PROGRAM 

STRUCTURE. THE RESULTS ARE COMPARED WITH DATA IN REF (C). 


DECLARATIONS... 
PARAMETER 

PARAMETER(MAX= 10000) 
ARRAYS 

REAL —A(3,3),L(3),UO(3),UI(3),U(3) 

DIMENSION TIME (MAX), ZBASE (MAX), FREQ (MAX), 
RELDIS (MAX), RELVEL (MAX), RELACC (MAX), 
EQUIPACC(MAX), EQMAXACC(MAX), ACCINPUT(MAX), 
ACCX1 (MAX), ACCX2 (MAX) 


OO.) Powree mm O00 OM OF © 


C VARIABLES 
REAL ZETA,DELF,DELT,DELT1,DELT2, TLOAD,E,F,G,W,MAXEQACC 
INTEGER NUMSTEP,NUMFREQ 

C PARAMETERS USED... 
ZETA 0. 


NUMFREQ = 250 
REWIND 10 
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2) CUS ekecec@, 


OO 


‘?) 


OO” "OCs ao 


OPEN(10,FILE=’TIMEREC .DAT,STATUS=’OLD’) 
OPEN(20,FILE=’"FREQREC .DAT,STATUS=’NEW’) 
READ AND LOAD THE BASE ACCELERATION ARRAY FROM THE DATA 


NOTE: (1) DELT EQUAL .0001 SECONDS 
(2) ACCELERATION IN G’S 
(3) FOR 2DOF TUNE FIXTURE, ACCX1 IS BASE ACCELERATION 
DO 100 I=1,NUMSTEP 
READ(10,*) TIME(I),ACCX1(I) 
READ(10,*) TIME(I),ACCX1(I),ACCX2(1) 
100 CONTINUE 
LOAD THE INPUT ACCELERATION ARRAY WITH THE BASE ACCELERATION IN 
G'S... 
DO 150 I=1,NUMSTEP 
ACCINPUT(I) =ACCX1(I)*32.2 
150 CONTINUE 
BEGIN ITERATIONS FOR THE NATURAL FREQUENCIES BEGINNING WITH 
DELF... 
DO 200 I=1,NUMFREQ 
FREQ(I)=I*DELF 


SOLVE E.O.M FOR INITIAL ACCELERATION WITH KNOWN PARAMETERS... 
UO(3) =-ACCINPUT(1) 
-(2.*ZETA*(2.*(4.*ATAN(1.))*FREQ(I)))*UO(2) 
-((2.*(4.* ATAN(1.))*FREQ(I))**2)  *UO(1) 


LOAD INITIAL VALUES FOR ITERATION AND COMPONENT ARRAYS.... 
Ul(1)  =UO(1) 7 
UI(2) =UO(2) 
UI(3)  =U0(3) 
RELDIS(1) =UO(1) 
RELVEL(1) =UO(2) 
RELACC(1) =U0(3) 


BEGIN TIME LOOP... 
DO 300 J=1,NUMSTEP 


SET TIME INTERVAL FOR PHASE 1 OR 2 OF TIME HISTORY....... 
IF((J*DELT1).LE.TLOAD)THEN 
DELT = DELT1 
ELSE 
DELT = 1./(20.*FREQ(I)) 
DELT2 = DELT 
ENDIF 


COMPUTE PARAMETERS FOR OPERATOR MATRICES... 
E=((2.*(4.*ATAN(1.))*FREQ(I))**2)*(DELT**2) 
F=(2.*ZETA*(2.*(4.*ATAN(1.))*FREQ(I)))*DELT 
G=(1./(1.+(F/2.)+(E/4.))) 


83 


SEe) 


Gre @ 


LOAD AMPLIFICATION AND LOAD VECTOR MATRIX... 
A(1,1)=G*(1.+(F/2.)) 
A(1,2)=G*DELT*(1.+(F/4.)) 
A(1,3)=G*DELT**2/4. 
A(2,1)=G*(-E/(2.*DELT)) 
A(2,2)=G*(1.-E/4.) 
A(2,3)=G*DELT/2. 
A(3,1)=G*(-E/((DELT**2)) 
A(3,2)=G*(-1.)*(E+F)/(DELT) 
A(3,3)=G*(-1.)*((F/2.)+(E/4.)) 
L(1) =G*(DELT**2)/4. 

L(2) =G*DELT/?2. 

L(3) =G 


BEGIN ITERATIONS... 
W =-ACCINPUT(J+1) 
U(1)= A(1,1)*UI(1)+A(1,2)*U1(2)+A(1,3)* U1(3)+W*L(1) 
U(2)= A(2,1)*UI(1)+A(2,2)* UI(2)+A(2,3)* U1(3)+W*L(2) 
UGB)= A@G,1)*UI(1)+AQG,2)*UI(2)+A(3,3)*UI13)+W*L(3) 


SAVE THE DIS,VEL,ACC ARRAY VALUES FOR THIS 


RELDIS(J+1) = U(1) 
RELVEL(J+1) = U(2) 
RELACC(J+1) = U(3) 


RESET THE ITERATION MATRIX FOR THE NEXT TIME 
ITERATION... 
DO 350 K=1,3 
UI(K)=U(K) 
350 CONTINUE 
300 CONTINUE 
END OF TIME ITERATION LOOP... 


COMPUTE THE EQUIPMENT ACCELERATION NOW THAT THE RELATIVE 
AND BASE ACCELERATION ARE KNOWN... 
DO 250 J=1,NUMSTEP 
EQUIPACC(J)=RELACC(J)+ACCINPUT(J) 
250 CONTINUE 


SELECT THE ABSOLUTE MAXIMUM FROM EQUIPACC ARRAY AND LOAD 
THAT VALUE INTO EQMAXACC ARRAY FOR THIS FREQUENCY ITERATE... 
CONVERT TO G’S... 
MAXEQACC=ABS(EQUIPACC(1)) 
DO 275 J=2,NUMSTEP 
IF(ABS(EQUIPACC(J)).GT.MAXEQACC)THEN 
MAXEQACC=ABS(EQUIPACC(J)) 
ENDIF 
275 CONTINUE 
EQMAXACC(I)=MAXEQACC#2.2 
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200 CONTINUE 
c END OF FREQUENCY LOOP... 


C UNITS FOR FREQUENCY RECORD ARE HZ AND G’S... 
DO 400 I=1,NUMFREQ 
WRITE(20,*) FREQ(I), EQMAXACC(I) 
400 CONTINUE 
END 
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APPENDIX D 


LT RANDALL CORBELL 

NAVAL POSTGRADUVATE SCHOOL 

SHOCK QUALIFICATION OF COMBAT SYSTEMS EQUIPMENTS USING TUNED 
MOUNTING FIXTURES ON THE U.S. NAVY MEDIUMWEIGHT SHOCK MACHINE 


PROGRAM: 2 DOF TUNED FIXTURE RESPONSE TO MWSM HALF-SINE 
ACCELERATION PULSE 


REF: (A) CLEMENTS,E.W.,"SHIPBOARD SHOCK AND NAVY DEVICES FOR 

ITS SIMULATION",NRL REPORT 7396, 1972. 

(B) CRAIG,R,"STRUCTURAL DYNAMICS, AN INTRODUCTION TO 
COMPUTER METHODS", JOHN WILEY AND SONS, 1981. 

(C) MEIROVITCH,L.,"“ELEMENTS OF VIBRATIONAL ANALYSIS", 
MCGRAW-HILL, 1986. 

(D) COSTANZO/MURRAY, "DTRC/UERD DDG-51 CLASS PRESHOCK 
TRIAL ANALYSIS PRELIMINARY REPORT’,1991. 

(E) COSTANZO/CLEMENTS,"MWSM CALIBRATION DATA", 1988. 


THIS PROGRAM COMPUTES THE RESPONSE OF A TWO DEGREE OF 
FREEDOM MASS-DAMPER-SPRING SYSTEM TO A BASE EXCITEMENTHALF-SINE 
ACCELERATION PULSE OF 1MSEC DURATION, AS CHARACTERISTIC WITH THE 
MWSM NOTED IN REF(A). REF(B) AND (C) WERE USED FOR NUMERICAL 
ANALYSIS AND PROGRAM STRUCTURE. THE UPPER TIER ACCELERATION 
RESPONSE OF THE 2 DOF TUNED FIXTURE IS THE BASE EXCITATION A 
COMBATSYSTEM EQUIPMENT WOULD EXPERIENCE IN UNDEX, AS SIMULATED 
ON THE MWSM. 


SELECTION OF THE COUPLED NATURAL FREQUENCIES OF THE SYSTEM ARE 
BASED ON ANALYSIS OF FFT/SHOCK SPECTRAL INFORMATION PROVIDED IN 
REF (D). SELECTION OF DAMPING VALUES WERE OBTAINED AFTER ANALYSIS 
OF MWSM CALIBRATION DATA, NOTED IN REF (BE). 


AS PER REF(B), A 2 DOF MASS SPRING SYSTEM SUBJECTED TO BASE 
MOTION CAN BE EXPRESSED IN TERMS OF RELATIVE MOTIONS: DISP, VEL, 
ACC AS FOLLOWS: 


DET Y1=X1-Z WHERE X1 IS MASS] COORDINATE AND Z IS BASE 
COORDINATE 
Y2=X2-Z WHERE X1 IS MASS2 COORDINATE AND Z IS BASE 
COORDINATE 
M1 = MASSI1 
M2 = MASS2 


K1 = SPRING STIFFNESS 1 
K2 = SPRING STIFFNESS 2 
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F1 = UNCOUPLED NATURAL FREQUENCY 1 
F2 = UNCOUPLED NATURAL FREQUENCY 2 
FN1= COUPLED NATURAL FREQUENCY 1 
FN2= COUPLED NATURAL FREQUENCY 2 
U11= MODAL MATRIX ELEMENT, (MASS1,FREQ1) 
U22= MODAL MATRIX ELEMENT, (MASS2,FREQ1) 
U12= MODAL MATRIX ELEMENT, (MASS1,FREQ2) 
U22= MODAL MATRIX ELEMENT, (MASS2,FREQ2) 
M1l= MASS MATRIX ELEMENT 
M21= MASS MATRIX ELEMENT 
M12= MASS MATRIX ELEMENT 
M22= MASS MATRIX ELEMENT 
K11= STIFFNESS MATRIX ELEMENT 
K21= STIFFNESS MATRIX ELEMENT 
K12= STIFFNESS MATRIX ELEMENT 
K22= STIFFNESS MATRIX ELEMENT 
Z11= BASE MATRIX ELEMENT 
Z21= BASE MATRIX ELEMENT 
Z12= BASE MATRIX ELEMENT 
Z22= BASE MATRIX ELEMENT 
DECLARATIONS... 
PARAMETER(MAX= 10000) 
ARRAYS... 
REAL _A(3,3),L(3), UO(3),UI(3), U(3) 
DIMENSION TIME (MAX),ANVILACC(MAX), 
ZBASE1 (MAX),ZBASE2 (MAX),ZBASE (MAX), 
RELDIS1(MAX),RELVEL1 (MAX), RELACC1(MAX), 
RELDIS2(MAX),RELVEL2 (MAX), RELACC2(MAX), 
YIACC (MAX),Y2ACC (MAX), 
: X1ACC (MAX),X2ACC (MAX) 
C VARIABLES... ) 
REAL F1,F2,FN1,FN2,W1,W2,WN1,WN2,W1SQ,W2SQ,WN1SQ,WN2SQ, 
M1,M2,M11,M21,M12,M22, 
K1,K2,K11,K21,K12,K22, 
U11,U21,U12,U22, 
Z11,Z21,ZCOEFF1,ZCOEFF2, 
ZETA,DELT,E,F,G,Z, WEIGHT 
INTEGER NUMSTEP,I1,12 
PRINT*, INPUT DAMPING RATIO, ZETA’ 
READ*, ZETA 
PRINT*, INPUT WEIGHT (LB) OF ONE TIER, MASS1=MASS2’ 
READ‘, WEIGHT 
PRINT*, INPUT NUMBER OF TIME STEPS UP TO 10,000” 
READ*, NUMSTEP 
C PARAMETERS USED... 


emer 1G) (ne) tC) () Crr) 0) 0 0 O 


Fl = 94, 

F2 = 100. 

FN1 = 60.691177 
FN2 = 154.88250 
DELT = 0.0001 
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vO(l) = 0. 

UO(2) = 0. 

M1 = WEIGHT/2.2 

M2 = M1 

Wl = 2.*4.*ATAN(1.)*F1 
w2 = 2.*4.*ATAN(I1.)*F2 
WNI1 = 2.*4.*ATAN(1.)*FN1 
WN2 = 2.*4.*ATAN(1.)*FN2 
W1SQ = W1**2 

Ww2SQ = W2**2 

WN1SQ = = WNI1**2 

WN2SQ. ss = WN2**2 

K1 = MI*W1SQ 

K2 = M2*W2SQ 


C NATURAL MODES... 
Ull = 1. 
U21 = (K1-WN1SQ*M1)/K1 
Ul12 = 1. 
U22 = (K1-WN2SQ*M1)K1 


Ss DECOUPLE SYSTEM... 
M11 = M1+U21*U21*M2 
M21 = M1+U22*U21*M2 
M12 = M1+U22*U21*M2 
M22 = M1+U22*U22*M2 


K11 = K1+U21*(-K1) + U21*(-K1+U21*(K1+K2)) 
K21 = K1+U21*(-K1) + U22*(-K1+U21*(K1+K2)) 
K12 = K1+U22*(-K1) + U21*(-K1+U22*(K1+K2)) 
K22 = K1+U22*(-K1) + U22*(-K1+U22*(K1+K2)) 


Z11 = M1+U21*M2 
Z21 = M1+U22*M2 


ZCOEFF1 = Z11/M11 
ZCOEFF2 = Z21/M22 


c READY TO PROCESS AS SDOF SYSTEMS... 


REWIND 10 
OPEN(10,FILE=’ANVILACCLM .DAT ,STATUS=’OLD’) 
OPEN(20,FILE=’TDOFTIMEREC.DAT ,STATUS=’NEW’) 


C READ DATA FROM TUNED FIXTURE FILE...ANVIL ACCELERATION IN G’S 
DO 10 I=1,NUMSTEP 
READ(10,*) TIME(I), ANVILACC(I) 
10 CONTINUE 


DO 20 I=1,NUMSTEP 
ZBASE (I) = ANVILACC(I)*32.2 
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ZBASE1(I) = ZCOEFF1*ANVILACC(I)*32.2 
ZBASE2(I) = ZCOEFF2*ANVILACC(I)*32.2 
20 CONTINUE 


C ITERATIONS... 
c 1. BEGIN WITH UPPER TIER SDOF MODEL...INITIALIZE PARAMETERS... 


UO(3) ee ee POR)WNISO"LIO() 


Ul(l) = 
Ul(2) = UO(2) 
UI(3) = UO(3) 


RELDIS1(1) = UO(1) 
RELVEL1(1) = UO(2) 
RELACC1(1) = UO(3) 


OPEN (30,FILE="TEST.DAT ,STATUS=’NEW’) 
l= 1 

WRITE(30,*) I1 

WRITE(30,*) F1,F2,FN1,FN2 

WRITE(30,*) W1SQ,W2SQ,WN1SQ,WN2SQ 
WRITE(30,*) M1,M2,K1,K2 

WRITE(30,*) U11,U12 

WRITE(30,*) U21,U22 

WRITE(30,*) M11,M12 

WRITE(30,*) M21,M22 

WRITE(30,*) K11,K12 

WRITE(30,*) K21,K22 

WRITE(30,*) Z11 

WRITE(30,*) Z21 

WRITE(30,*) K11/M11,K22/M22 

WRITE(30,*) ZCOEFF1 

WRITE(30,*) ZCOEFF2 

WRITE(30,*) UO(1),UO(2),UO0(3) 

WRITE(30,*) RELDIS1(1), RELVEL1(1),RELACC1(1) 


cS COMPUTE PARAMETERS FOR OPERATOR MATRICES... 
E=WN1SQ*(DELT**2) 
=2.*ZETA*WN1*DELT 
G=(1./(1.+(F/2.)+(E/4.))) 


c LOAD AMPLIFICATION AND LOAD VECTOR MATRIX... 
A(1,1)=G*(1.+(F/2.)) 
A(1,2)=G* DELT* (1.+(F/4.)) 
A(1,3)=G*DELT**2/4. 
A(2,1)=G*(-E/(2.*DELT)) 
A(2,2)=G*(1.-E/4.) 
A(2,3)=G*DELT/2. 
A(3,1)=G*(-E/((DELT**2)) 
A(3,2)=G*(-1.)*(E+F)(DELT) 
A(3,3)=G*(-1.)*((F/2.) +(E/4.)) 
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L(1) =G*(DELT**2)/4. 
L(2) =G*DELT/2. 
L3) =G 


BEGIN TIME ITERATIONS UPPER TIER MOTION... 
DO 100 J=1,NUMSTEP 

Z =-ZBASE1(J+1) 

U(1)= A(1,1)*UI(1)+A(1,2)*UI(2)+A(1,3)*UI(3)+Z*L(1) 
U(2)= A(2,1)*UI(1)+A(2,2)*UI(2)+A(2,3)*UI(3)+Z*L (2) 
U(3)= A(3,1)*UI(1)+A(,2)*UI(2)+AG,3)*UI(3)+Z*L(3) 


SAVE THE DIS,VEL,ACC ARRAY VALUES FOR THIS ITERATION... 
RELDIS1(J+1) = U(1) 
RELVEL1(J+1) = U(2) 
RELACC1(J+1) = U(3) 


RESET THE ITERATION MATRIX FOR THE NEXT TIME 
ITERATION... 

UI(1)=U(1) 

UI(2)=U(2) 

UI(3)=U(3) 

100 CONTINUE 


2. NOW FOR LOWER TIER SDOF MODEL...INITIALIZE PARAMETERS... 


UO(3) =-ZBASE2(1)-(2.*ZETA*WN2)*UO(2)-WN2SQ*UO(1) 
Ul(1) = UO(1) 
Ul(2) = UO(2) 
UI(3) = UO(3) 


RELDIS2(1) = UO(1) 
RELVEL2(1) = UO(2) 
RELACC2(1) = UO(3) 

I2=2 

WRITE(30,*) 12 

WRITE(30,*) F1,F2,FN1,FN2 
WRITE(30,*) W1SQ,W2SQ,WN1SQ,WN2SQ 
WRITE(30,*) M1,M2,K1,K2 
WRITE(30,*) U11,U12 
WRITE(30,*) U21,U22 
WRITE(30,*) M11,M12 
WRITE(30,*) M21,M22 
WRITE(30,*) K11,K12 
WRITE(30,*) K21,K22 
WRITE(30,*) Z11 

WRITE(30,*) Z21 

WRITE(30,*) K11/M11,K22/M22 
WRITE(30,*) ZCOEFF1 
WRITE(30,*) ZCOEFF2 
WRITE(30,*) UO(1),UO(2),UO(3) 
WRITE(30,*) RELDIS2(1), RELVEL2(1),RELACC2(1) 


Ooo 


200 


COMPUTE PARAMETERS FOR OPERATOR MATRICES... 
E=WN2SQ*(DELT**2) 

F=2.*ZETA*WN2*DELT 

G=(1./(1.+(F/2.)+(E/4.))) 


LOAD AMPLIFICATION AND LOAD VECTOR MATRIX... 
A(1,1)=G*(1.+(F/2.)) 
A(1,2)=G*DELT*(1.+(F/4.)) 
A(1,3)=G*DELT**2/4. 
A(2,1)=G*(-E/(2.*DELT)) 
A(2,2)=G*(1.-E/4.) 
A(2,3)=G*DELT/?2. 
A(3,1)=G*(-E/(DELT**2)) 
A(3,2)=G*(-1.)*(E+F)/(DELT) 
A(3,3)=G*(-1.)*((F/2.)+(E/4.)) 
L(1) =G*(DELT**2)/4. 

L(2) =G*DELT/2. 

L(3) =G 


BEGIN ITERATIONS... 

DO 200 J=1,NUMSTEP 

Z =-ZBASE2(J+1) 

U(1)= A(1,1)*UI(1)+A(1,2)*UI(2)+A(1,3)*U1(3)+Z*L(1) 
U(2)= A(2,1)*UI(1) +A(2,2)*UI(2) +A(2,3)*UI(3)+Z*L (2) 
U(3)= A(3,1)*UI(1) +A(3,2)*U1(2) +A(3,3)*UI(3)+Z*L(3) 


SAVE THE DIS, VEL,ACC ARRAY VALUES FOR THIS ITERATION... 
RELDIS2(J+1) = U(1) 
RELVEL2(J+1) = U(2) 
RELACC2(J+1) = U(3) 


RESET THE ITERATION MATRIX FOR THE NEXT TIME 
ITERATION... 

UI(1)=U(1) 

UI(2)=U(2) 

UI(3)=U(3) 

CONTINUE 


CONVERT RELATIVE UNCOUPLED MOTIONS TO RELATIVE COUPLED MOTIONS 
USING MODAL MATRIX...THEN RELATIVE COUPLED MOTIONS TO ABSOLUTE 
TIER MOTIONS... 


300 


DO 300 I=1,NUMSTEP 
YIACC(I)= U11*RELACCI(I) + U12*RELACC2(I) 
Y2ACC(I)= U21*RELACCI(I) + U22*RELACC2(I) 
XIACC(I)= (YIACC(I) —- + ZBASE(I))2.2 
X2ACC(I)= (Y2ACC(I) + ZBASE(I))32.2 
WRITE(20,*) TIME(I),X1ACC(I), X2ACC(I) 

CONTINUE 

END 


91 


Oe GG Gre @ Ge Goma. 70 OO -0OO 


2) ie 


APPENDIX E 


LT RANDALL CORBELL 
NAVAL POSTGRADUATE SCHOOL 


SHOCK QUALIFICATION OF COMBAT SYSTEMS EQUIPMENTS USING TUNED 
MOUNTING FIXTURES ON THE U.S. NAVY MEDIUMWEIGHT SHOCK MACHINE 


2 DOF TUNED FIXTURE COUPLED NATURAL FREQUENCIES STUDY FOR 
MASS1/MASS2 RATIO 


REF: (A) SHIN, Y.S.,"NEDE REPORT NED345, CLASS IT", 1981 


THIS PROGRAM ITERATES THE EXPRESSIONS FOR THE COUPLED NATURAL 
FREQUENCIES, FN1 AND FN2, OF A 2 DOF SYSTEM AS EXPRESSED IN TERMS 
OF ITS UNCOUPLED NATURAL FREQUENCIES, F1,F2, AND THE MASS RATIO 
OF THE SYSTEM, MASS1/MASS2. REF(A) WAS USED TO OBTAIN AN 
EXPRESSION FOR THIS RELATION. 

FN1 = FIRST MODE COUPLED NATURAL FREQUENCY 

FN2 = SECOND MODE COUPLED NATURAL FREQUENCY 

F12 = FIRST MODE UNCOUPLED NATURAL FREQUENCY SQUARED 

F22 = SECOND MODE UNCOUPLED NATURAL FREQUENCY SQUARED 

R = MASS1/MASS2 


DECLARATIONS... 
REAL F1,F2,F12,F22,FN1,FN2,R,DELF,C 
PRINT*, INPUT MASS1/MASS2 RATIO, R’ 
READ’, R 
PRINT*, INPUT FI’ 
READ*, Fl 
DELF=1. 
C =(5)**.5 
OPEN(20,FILE=’NATFREQDAT.DAT,,STATUS=’NEW’) 
DO 100 I=1,250 
Fl = I*DELF 
Fi2 = Fi**2 
DO 200 J=1,400 
F2 = J*DELF 
F2 = 100. 
F22 = F2**2 
FN1 = C*((F12+R*F12+F22) 
-((F12+ R*F12+F22)**2 - 4.*F12*F22)**.5)**.5 
FN2 = C*((F12+R*F12+F22) 
+ ((F12+R*F12+F22)**2 - 4.*F12*F22)**.5)**.5 
IF(((FN1.GE. 20.).AND.(FN1.LE. 25.)).AND. 
((FN2.GE.153.).AND.(FN2.LE.156.))) THEN 
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re) 


WRITE(20,*) F1,F2,FN1,FN2,R 


ENDIF 
200 CONTINUE 
100 CONTINUE 
END 
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APPENDIX F 


The preliminary design for a Two Degree of Freedom MWSM Tuned 
Mounting Fixture is presented below. The beam loading and configuration 
determines the system stiffness. For this model, each tier is simply supported with 
uniform loading over the span. L, between the two simple support beams of each 


tier. Figure F-1 depicts the model and Table F-1 lists the parameters. 


TABLE F-1 


PARAMETER 
NATURAL FREQUENCY HZ ae a 

| WEIGHT LBS 
MASS LB*SEC~2/FT 
STIFFNESS LB/FT 
BEAM LENGTH W5x18.5 FT 
NO. STIFFNESS BEAMS 
NO. SIMPLE SUPPORT BEAMS 
SUPPORT BEAM SPACING,LFT | 40 | 38 — 


| LOWER TIER MOUNTING PLATE 12 
THICKNESS IN 



















il 






94 


2 DOF TUNED MOUNTING FIXTURE PRELIMINARY DESIGN 


SPECIFICATIONS: UPPER TIER 1700 LBS, F= 94HZ 
LOWER TIER 1700 LBS, F = 100HZ 
BEAMS, AISC W5 x 18.5 





UPPER TIER 
LOWER TIER 
fa ee 
512 ae 1 ZL, 
A me 
ANVIL TABLE 5XSFT 
Figure F-1. Preliminary Design for Two Degree of Freedom Tuned 


Mounting Fixture. 
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